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A model which describes the time-dependent CP formalism in D° decays has recently been 
proposed. There it has been highlighted a possible measurement of the angle f3 c , in the charm 
unitarity triangle, using the decays D° -> K + K~ and D° — > 7r + 7r , and a measurement of the 
mixing phase 4>mix- The same method can be used to measure the value of the parameter x, one of 
the two parameters defining charm mixing. We numerically evaluate the impact of a time-dependent 
analysis in terms of the possible outcomes from present and future experiments. We consider the 
scenarios of correlated D° mesons production at the center of mass energy of the ^(3770) at SuperB, 
uncorrelated production at the center of mass energy of the T(4S) at SuperB and Belle II, and LHCb. 
Recently a hint of direct CP violation in charm decays was reported by the LHCb collaboration, 
we estimate the rate of time-dependent asymmetry that could be achieved using their available 
data, and we generalise the result for the full LHCb program. We conclude that LHCb is already 
able to perform a first measurement of /3 c ,e//, and slightly improve the present constraints on the 
parameters x and 4>mix- A more precise determination of /3 c ,e//, ^mix and x will require a larger 
data sample, and most probably the cleaner environment of the new high luminosities B-factories 
(both SuperB and Belle II) will be needed. We show that SuperB will be able to measure /9 c ,e// 
and 4>mix with a precision of 1.4° and improve the precision on i by a factor of two. 

PACS numbers: 13.25.Hw, 12.15.Hh, 11.30. Er 



I. INTRODUCTION 

Since the discovery in 1964 of CP violation in the Kaon 
system li , CP violation has been observed also in the B 
meson system [2] [3]. In the charm sector, CP violation 
has long been expected to be too small to be observed 
at precision available until recently when, the LHCb col- 
laboration has reported a difference in direct CP asym- 
metries in D° K + K~ and D° -> 7r+7r~ that is 3.5ct 
from the CP conserving hypothesis [I]. In [5] the stan- 
dard model (SM) description of these decays using the 
same CKM paradigm that provides a rather satisfactory 
description of such decays of B° mesons is considered. 
Since the LHCb result, a broader view of this paradigm 
that might accomodate the large asymmetry is examined 
in [6]. It is clear that, in order to understand the nature 
of CPVin D° decays, measurements of weak phases in 
these decays are essential. In [5], it is proposed that, as 
with B° decays, time-dependent CP asymmetries in D° 
decays may provide the most direct way to measure these 
phases. In this paper, we further examine the precision 
that might be anticipated in four experimental scenarios 
that are likely to be available over the coming decade to 
evaluate the rate of time-dependent CP asymmetries in 
D° — > K + K~ and D° — > ir + ir~ in the three proposed 
environments (SuperB, LHCb, Belle II). 



II. TIME-DEPENDENT CP VIOLATION IN 
THE CHARM SECTOR 

In the standard model (SM), CP violation is described 
in terms of the complex phase appearing in Cabibbo- 
Kobayashi-Maskawa (CKM) matrix [7] [5]. The matrix 
is a unitary 3x3 matrix which provides a description of 
quark mixing in terms of the coupling strengths for up- 
to-down quark type transitions, and it may be written 
as 

(V uc i V us V u b \ 
v cd V cs V cb . (1) 
V td V u V tb J 

Within this framework the probability to observe a tran- 
sition between a quark q to a quark q' is proportional to 
\V qq >\ 2 - 



A. Buras parametrisations of the CKM Matrix 

In Ref. [5] the CKM matrix has been written using the 
Buras" parametrisation 



l-A 2 /2-A 4 /8 A AX 3 (p-ifj) + AX 5 {p-ifj)/2 

Vckm = ( -A + A 2 A 5 [1-2Q5 + if})}/2 1 - A 2 /2 - A 4 (l + 4A 2 )/8 AX 2 |+0(A H ),(2) 

AA 3 [1 - (p + ifj)] -AX 2 + AX i [l-2(p + ifj)}/2 l-A 2 X i /2 
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We adopt the convention of writing the CKM matrix 
in terms of p and rj because these represent the coordi- 
nates of the apex of the well known bd unitarity triangle. 
Since unitarity triangles are mathematically exact, it is 
very important to measure their angles and sides to ver- 
ify unitarity. One of the six unitarity relationships of the 
CKM matrix may be written as 

v: d v cd + v: s v cs + v; b v cb = o, (3) 

where Eq. ([3| represents the cu triangle that we will call 
the charm unitarity triangle or simply charm triangle. 
The internal angles of this triangle are given by 

a c = Mg[-V: b V cb /V: s V cs ], (4) 

& = Mg[-v; d v cd /v: s v cs ], (5) 

7c = Mg[-V; b V cb /V: d V cd }. (6) 

In Ref. [5] we proposed the measurement of fi c ,eff using 
time-dependent CP asymmetries in charm decays and 



using the results of Global CKM fits, predicted that: 

(3 C = (0.0350 ±0.0001)°. (7) 

On comparing Eq. ^ with Eq. one can see that 
V cd = V c de l ^ a ~^ in this convention. 



B. Time-dependent formalism 

We consider two different cases of D° meson pro- 
duction: un-correlated and correlated D° production. 
Un-correlated D 0, s are produced from the decays of B 
mesons in electron-positron colliders when particles are 
collided at a center of mass energy corresponding to the 
T(AS) resonance, from cc continuum, or in hadrons col- 
lider. The correlated D° mesons are produced in an 
electron-positron machine running at a center of mass en- 
ergy corresponding to the '5 (3770) resonance. The time 
evolution for both situations is given by [5] 



un-correlated case 



T+ oc e 



-Tit 



1 + e 
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l + IA/l 2 
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(8) 



correlated 



T + cx e- ri l A *l 



h+ Re(Xf) 
2 + l + \\ f \* 



h _ ± e ArAt/2 



l + IA/l 5 
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2Im(\ f ) 
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sin AM At 



(9) 




FIG. 1: Generated distributions according to our formulas for 
D° — > /(left) and for D° — > /(right) produced at the center- 
of-mass energy of the $(3770). 



to 



Where T + refers to D°(q c = +2/3) decays and T_ 
W{q c = -2/3) decays, h± = l± e ArAt and A/ = 
Here q and p are the parameters defining the mixing and 



A (A) is the amplitude for the D (D) decay to a final 
state /. If \A\ 2 ^ \A\ 2 there is direct CP violation (in 
the decay) and \q/p\ =/= 1 would signify CP violation in 
mixing. The study of A/ (which should not be confused 
with the term A appearing in the CKM matrix) is able to 
probe the combination of CP violation due to mixing and 
decay, and this form of CP violation is referred to as CP 
violation in the interference between mixing and decay. 
Considering Eqns. (8j9) the time-dependent asymmetries 
associated with the time evolution of the D° mesons can 
be written in terms of the physical decay rate including 
the mistag probability, u}(u>), for incorrect tagging of the 
D° (D°) flavour as follows 



T Phys {t) = (i_ w)r+ (t) + ajr_(<), 

cjr_i 



T Phys (t) 



(10) 

+ (f) + (i-w)r_(t), (ii) 

where r+(i) and r_(£) are from Eqns. ^ and ^ and 
uj (to) represents the mistag probability for the particle 
(anti-particle) apparent decay rates for D° and D°, re- 
spectively. Hence for un-correlated mesons the time de- 
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pendent CP asymmetry accounting for mistag probabil- ity is 



A Phys {t) = 



-Phys 



(t)-T ph v s (t) 

T Phvs {t) + r Ph v 3 (t) 



= Au 



(D - Aw)e Art / 2 [(|A/| 2 - l)cosAM£ + 2Im\ f sinAMt] 
h + (l + \\ f \)y2 + Re(\ f )h- : 



(12) 



where Aw 



oj and D = 1 — 2u>. 



Similarly the asymmetry for correlated mesons is 



A Phys (At) = 



^Phys 
-Phys 



r '""(At) - r Ph y s (M) 



= -Aw ■ 



r " (At) + T Ph v*(At) 



(D + Aw)e ArAt / 2 [(|A / | 2 - 1) cos AAfA< + 27mA/ sinAMAt] 
h + (l + \X f \)y2 + Re(X f )h^ 



(13) 



The above equations may be written in terms of x and 
y allowing for the measurement of the mixing phase. We 
report here the time- dependent asymmetry equation for 



correlated mesons (similar results may be obtained in the 
un-correlated case): 



A Ph y y s (At) = -Au 



(D + Aw)e y [(|A/| 2 - l)cosa;rAt + 2/TOA / sinzTAt] 
h+(l + \\ f \)y2 + Re(\ f )h_ • 

I 



(14) 



III. MC TEST OF TIME-DEPENDENT 
NUMERICAL ANALYSIS 



One of the issues raised in [5] is the possibility to use 
different decay channels of the D° mesons to constrain 
the value of the angle f3 c of the charm triangle. The de- 
cay D° — > K + K~ will be used to measure the mixing 
phase, the decay D° — > 7r + 7r~ will be used to measure 
<Pmix — 2/3 c and the difference between the two channels 
will provide a first measurement of the angle j3 c . In this 
framework, long distance contributions to decay are not 
considered. The latter together with the different con- 
tribution to decay D° — > tt + tt^ from penguin topologies 
will introduce theoretical uncertainties, and for this rea- 
son we refer to the angle (3 C as Pc,eff where effective in- 
dicates that there are theoretical uncertainties that need 
to be evaluated. To evaluate the asymmetry, and es- 
timate the precision on c ,eff that one might achieve 
in the different experimental environments described in 
the previous section, we generate a set of one hundred 
Monte Carlo data samples. Each one based on the ex- 
pected number of tagged D° decays in the corresponding 
experimental setup, and we generate data according to 
the distributions given in Eqns. ^ and ([9|, where the 
parameters involved are evaluated as in Ref. 11 . We 



evaluate the asymmetry given in Eqns. (12 1 and (13) in 



eluding the expected mistag probabilities, and perform a 
binned fit to the simulated data. The distributions that 
we are considering here have been expressed as functions 
of | Ay | and arg(Xf) = <j> = 4>mix — 2<PcPi and the fit 
is performed keeping |A/| — 1 and allowing arg(Xf) to 
vary. The same results are obtained when also |A/| is 
also allowed to vary in the fit. It is important to mention 
that a measurement of Ay 7^ 1 in an experiment would 
be a signature of CP violation [5]- 



A. SuperB at the T(4S) 

The Superi? collaboration is expected to start taking 
data in 2017 PJ] [B] [H] [TS], and the integrated lu- 
minosity which will be achieved with the full program is 
expected to be 75 ab -1 . With this luminosity one would 
expect to reconstruct 6.6 x 10 6 tagged D° —> n + ir~ events 
in a data sample of 75 ab -1 with a purity of 98% [S]. The 
results of the numerical analysis are shown in Fig. [2] 

The asymmetry parameters determined here have a 
precision of (Jargix,,,,) — = 2.2°. The same pro- 

cedure when applied to the D° —> K + K~ channel to 
measure <J a rg(x KK ) = a 4>KK^ f° r which one would expect 
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FIG. 2: The time-dependent CP asymmetry expected for 
D° — > it + tt~ decays in a 75 ab _1 sample of data at the T(4S). 



to reconstruct 1.8 x 10 7 events, leads to precision of 
a ( f >KV — 1.6°. When the results from D° — > K + K~ and 



D° 



are combined one obtains a precision in 
, = 1.4°. 



B. SuperB at the *(3770) 

The Superi? collaboration is planning to have a ded- 
icated run at the center of mass energy of the \& (3770) 
resonance, to collect an integrated luminosity of 1.0 ab _1 . 
With this luminosity one would expect to record 979000 
D° —> 7r + 7r~ reconstructed events, when the full set of 
semi-leptonic decays K^ivg I = e, [i is used to tag the 
flavor of D° mesons (with negligible mistag probability). 
The results of the numerical analysis are shown in Fig. [3] 



X is anything, which corresponds to 54% (3%) of all D° 
meson decays from which one would expect to ~ 0.03, 
and that the asymmetry in particle identification of K + 
and K~ in the detector will naturally lead to a small, 
but non-zero value of Aw. We expect that there would 
be approximately 4.8 million kaon tagged D° — > ir + ir~ 
events in l.Oab -1 at charm threshold. Using these data 
alone, one would be able to measure cf) n7r to a precision 
of 2.7°. Hence if one combines the results from semi- 
leptonic and kaon tagged events, a precision of u^^^ ~ 
2.4° is achievable. 



C. LHCb 

Another possible scenario is that of measuring time- 
dependent asymmetries from uncorrelated D mesons in 
a hadronic environment, in particular the LHCb experi- 
ment. Here dilution and background effects will be larger 
than those at an e + e~ machine, but the data are al- 
ready available and it would be interesting to perform 
the time-dependent analysis, especially after the recent 
results on time integrated CP violation in Ref. [J]. As 
already mentioned a measurement of |A/| ^ 1 will signify 
CP violation. Given that the measurement of A/ is likely 
expected to be dominated by uncertainties, especially in 
lj and Aw, it is not clear what the ultimate precision 
obtained from LHCb will be. The best way to ascertain 
this would be to perform the measurement on the exist- 
ing data set. We have estimated that LHCb will collect 
4.9 x 10 6 D* tagged D° ->■ tt+tt" decays in 5fb _1 of 
data, based on the 0.62 fb" 1 of data shown in [4], and we 
consider also the outcome of a measurement for l.lfb -1 
(equivalent to 0.7 x 10 6 D* tagged D° -» ir + n~ decays) 
already available after the 2011 LHC run. In [5] we esti- 
mate a purity of ~ 90% and lu ~ 6% which results in the 
asymmetry obtained in Fig. [3] for 5 fb _1 of data. 




FIG. 3: The time-dependent CP asymmetry expected for 
D° — > 7r + 7r~ decays in 1 ab _1 sample of data at the ^(3770). 
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FIG. 4: The time-dependent CP asymmetry expected for 
D° — > 7r + 7r~ decays in a 5fb _1 sample of data at LCH6. 



The phase 0^ could be measured with a precision of 
= 5.7°. One may also consider using hadronically 
tagged events, for example D° —> K~X (K + X), where 



This fit is translated into a potential measurement of 
the phase with a precision of 3.0°. With l.lfb -1 
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of data we estimate that LHCb may be able to reach a 
precision of 8° on 4>-mr- 

D. Belle II 

The last scenario considered here is that of Belle II with 
50 ab -1 of data collected at the center of mass energy of 
the Y(AS) [IB]. We have considered the same efficiency 
and mistag probability as for SuperB and we expect that 
4.4 x 10 6 D* tagged D° -> tt+tt" will be collected. The 
resulting asymmetry is shown in Fig. [5} The precision on 
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FIG. 5: The time-dependent CP asymmetry expected for 
D° — > tt+tt~ decays in a 50 ab -1 sample of data at the T(4S) 
at Belle II. 



obtained for this scenario is estimated to be 2.8°. 

IV. TIME-DEPENDENT SENSITIVITY 
STUDIES 

A. Sensitivity to x 

We consider the same data sample discussed in the pre- 
vious sections for D° —> tt+tj- and D° -+ K+K- . While 
we find that results from the time-dependent analysis are 
not sensitive to the parameter y, and that with 1.0 ab -1 
of data collected at charm threshold at SuperS it will be 
possible to improve the currently known precision on x 
by a factor of two with respect to the most recent HFAG 
values [17] . The precision that could be reached is shown 
in Table E 

B. Sensitivity to Pc,eff, </>mix and (/>cp 

We show here a summary of the possible sensitivities 
that the different experiments could achieve when mea- 
suring the mixing and the weak phase. 

At first order the decay D° —> K + K~ measure 
the mixing phase, therefore one can consider <f>KK = 



TABLE I: Estimates of the sensitivity on x for all the exper- 
imental scenarios and their projected luminosities for the de- 
cays D° — > tt+ty~ and D° — » K+K~ and <j> = cj)Mlx — 2Pc,eff- 



Experiment/HFAG 


a w (4> = ±10°) 


(7,(0 = ±20°) 


SuperB [T(45)] 






D° -> 7T+7T- 


0.12% 


0.06% 


D° -> K+K" 


0.08% 


0.04% 


SuperB [#(3770)] 






D° -> 7r+7r"(S'L) 


0.30% 


0.15% 


D° -» n+n (SL + K) 


0.13% 


0.06% 


D° -> K+K-(SL) 


0.19% 


0.10% 


D° -» K+K" (SL + K) 


0.08% 


0.04% 


LHCb 






D° -» 7r + 7r (1.1 fb" 1 ) 


0.40% 


0.20% 


D° -t K+K- (1.1 fb _1 ) 


0.22% 


0.11% 


D° -» TT+TY- (5.0 ftT 1 ) 


0.15% 


0.08% 


D° -> K+K- (5.0 fb _1 ) 


0.09% 


0.04% 


Belle II 






D° -> ir+TT- 


0.14% 


0.07% 


D° -> K+K~ 


0.10% 


0.04% 


HFAG 


0.20% 



TABLE II: Summary of expected uncertainties from 1 ab 
of data at charm threshold, 75 ab -1 of data at the T(4S), 
5fb _1 of data from LHCb, and 50 ab -1 of data at the Y(4S) 
at Belle II for D° —¥ tt+it~ decays. The column marked SL 
corresponds to semi-leptonic tagged events, and the column 
SL+K corresponds to semi-leptonic and kaon tagged events at 
charm threshold. The last row shows the precision in /3 c ,e// 
expected from a simultaneous fit to mr and KK where we 
assume that, for KK, the decay is dominated by a tree am- 
plitude. 







SuperB 


LHCb 


Belle II 


Parameter 


#(3770) #(3770) T(4S) 


nf 






SL 


SL+K 7v+ 


nt 


= cr arg(X„„) 


5.7° 


2.4° 2.2° 


3.0° 


2.8° 


a <t>KK a arg(X K K ) 


3.5° 


1.4° 1.6° 


1.8° 


1.8° 




3.3° 


1.4° 1.4° 


1.9° 


1.7° 



arg(\KK) — 4>mix and use the time dependent analysis 
to measure it to a precision of 1.4° — 1.6°. 



C. Systematic uncertainties 

The knowledge of the parameters x and y which de- 
fine the mixing is limited by their relative uncertainties. 
Since our analysis is not sensitive to the parameter y, we 
considered the most recent results from the HFAG [17] 
and we evaluated the effect of varying the parameter 
Ar = 2yT considering plus-and-minus one standard de- 
viation. This is the systematic uncertainty due to the 
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limited precision in y. The value of the uncertainty in 
the parameter y is 0.013% and it is given in [17] . The 
results are shown in Table IIIII 

TABLE III: Summary of expected systematic uncertainties 
due to the limited knowledge of the parameter y from 1 ab~ 
of data at charm threshold and 75 ab _1 of data at the T(4S). 
The column marked SL corresponds to semi-leptonic tagged 
events, and the column SL+K corresponds to semi-leptonic 
and kaon tagged events at charm threshold while nf 1 refers to 
the slow pion tag at the T(4S I ). 



Parameter 


*(3770) *(3770) T(4S) 
SL SL+K Tvf 


°<t>K K { s y s -) 
<7ts c , eIf {sys.) 


0.5° 0.2° 0.05° 
0.2° 0.1° 0.02° 
0.27° 0.11° 0.03° 



D. Combined results for SuperB 

We evaluated the combination of the results obtained 
for the different centre of mass energy at SuperB. The 
final results are made on the assumption that cf> = 
4>mix — 2/3 c = ±10° and they are shown in Table |TV| 



V. CONCLUSIONS 



This paper elucidates the time-dependent analysis of 
the D° mesons discussed in Ref. [5]. We concentrated 
on the possible measurement of the (3 c .eff angle of the 
charm unitarity triangle, on the mixing phase (j)Mix and 
on the mixing parameters. We estimate our results and 
compare them for the experimental environments that we 
think could and should perform this analysis: SuperB, 
LHCb and Belle II. We found that SuperB may perform 
better this analysis, but time is required before the col- 
laboration will start data taking. LHCb will have to con- 
trol the background levels to perform this measurement 
resulting then in a challenging analysis. However as re- 
ferred to in the article the LHCb collaboration has al- 
ready available an amount of data to analyse. This same 
amount of data has already shown a first hint of direct 
CP violation in charm, we think it would be worth go 
through the time-dependent formalism. The Belle II col- 
laboration will start data taking in few years, and the 
background-clean environment will allow to perform a 
time-dependent analysis and an evaluation of the mixing 
phase and of the /3 c ,e// at high precision. 



TABLE IV: Combined sensitivities at SuperB. 



Parameter 


Statistical 


Systematic 


sensitivity 


sensitivity 


a x (B u ->■ 7T+7T-) 


0.09% 




a x (D° -> K+K-) 


0.05% 




^„ 


1.62° 


0.14° 


a <l>KK 


1.05° 


0.02° 




0.92° 


0.03° 
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